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ABSTRACT 

We present a deep Chandra observation of the central regions of the late-type barred spiral galaxy 
NGC 2903. The Chandra data reveal soft (kT e ~ 0.2-0.5 keV) diffuse emission in the nuclear starburst 
region and extending ~2' (~5 kpc) to the north and west of the nucleus. Much of this soft hot gas is 
likely to be from local active star-forming regions; however, besides the nuclear region, the morphology 
of hot gas does not strongly correlate with the bar or other known sites of active star formation. The 
central ^650 pc radius starburst zone exhibits much higher surface brightness diffuse emission than the 
surrounding regions and a harder spectral component in addition to a soft component similar to the 
surrounding zones. We interpret the hard component as being also of thermal origin with kT e ~3.6 keV 
and to be directly associated with a wind fluid produced by supernovae and massive star winds similar 
to the hard diffuse emission seen in the starburst galaxy M82. The inferred terminal velocity for this 
hard component, ^1100 km s _1 , exceeds the local galaxy escape velocity suggesting a potential outflow 
into the halo and possibly escape from the galaxy gravitational potential. Morphologically, the softer 
extended emission from nearby regions does not display an obvious outflow geometry. However, the 
column density through which the X-rays are transmitted is lower in the zone to the west of the 
nucleus compared to that from the east and the surface brightness is relatively higher suggesting some 
of the soft hot gas originates from above the disk; viewed directly from the western zone but through 
the intervening disk of the host galaxy along sightlines from the eastern zone. There are several 
point-like sources embedded in the strong diffuse nuclear emission zone. Their X-ray spectra show 
them to likely be compact binaries. None of these detected point sources are coincident with the mass 
center of the galaxy and we place an upper limit luminosity from any point-like nuclear source to be 
<2xl0 38 erg s _1 in the 0.5—8.0 keV band which indicates that NGC 2903 lacks an active galactic 
nucleus. Heating from the nuclear starburst and a galactic wind may be responsible for preventing 
cold gas from accreting onto the galactic center. 

Subject headings: galaxies: evolution — galaxies: individual (NGC 2903), — galaxies: nuclei - 
X-rays: galaxies 



1. INTRODUCTION 

Stellar bars are common features of disk galaxies (Mas- 
ters et al. 2011; Sellwood 1993; Buta 2011) where they 
efficiently channel gas into the central regions. The dy- 
namics of bars is critical to galactic (pseudo)bulge growth 
and hence to evolution along the Hubble sequence (Kor- 
mendy & Kennicutt 2004) and to fueling central super- 
massive black holes (e.g., Shlosman et al. 1990) in iso- 
lated galaxies. In isolated late-type disk galaxies, stellar 
bars can efficiently transport gas into their central re- 
gions because of their relatively low bulge masses com- 
pared to early-type disks. Gas may accumulate at the 
ends of the bar before being driven inward by the grav- 
itational torques within the bar (Combes & Elmegrecn 
1993). The gas forms thin "centered" dust lanes along 
the bar (Athanassoula 1992). There is less shear in the 
bars of late-type disks giving rise to substantial star for- 
mation along the bar. 

There are often no central resonances in late-type 
disks, in contrast to bars in early-type disks; the bar ex- 
tends to about the turn-over radius of the rotation curve, 
near the inner Lindblad resonance, and corotation oc- 
curs even further out in the disk (e.g., Shlosman 1999). 
In this case, accretion does not stall to form a circum- 



nuclear ring but accumulates until self-gravity leads to 
burst-like episodes (Sarzi et al. 2007) of star formation - 
inflow may be impeded by stellar feedback but not by 
resonance phenomena as in early-type disks. Star forma- 
tion and feedback can deplete the gas reservoir even at 
small distances beyond the sphere of influence of central 
black holes (Davies et al. 2010; Thompson et al. 2005) or, 
it may enhance fueling of black holes if the starburst oc- 
curs very near the central compact object (Watabe et al. 
2008). 

The study of structures such as bars and circumnu- 
clear rings in X-ray light can directly reveal details that 
are hidden at many other wavelengths. The X-ray mor- 
phology not only helps to pinpoint locations of recent 
star formation but also traces high-pressure hot gas wind 
fluid that may quench star formation, impede accretion 
toward the galactic center, and drive energetic outflows. 
The best subjects for X-ray studies of these phenomena 
are isolated nearby galaxies inclined enough to distin- 
guish galactic outflows while leaving the nucleus exposed 
to direct view and galaxies lacking strong active galac- 
tic nucleus (AGN) activity which could mask the fainter 
X-ray emission from circumnuclear star formation. 

NGC 2903 is a nearby (8.9 Mpc, 1"=43 pc; Drozdovsky 
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Fig. 1. — Left: the central 4'x 6' region of SDSS image. Blue, green, and red correspond to the g, r, and i bands, respectively. Young 
stars in spiral arms are traced in blue, and bright Hn regions appear in green. The complex dust lane structures along the bar and spiral 
arms are also noticeable in darker color. White box corresponds to the region shown in the right panel. Right: the central 2'x3' region of 
Chandra image. Red, green, and blue correspond to the energy ranges, 0.5—1.0 keV, 1.0-2.0 keV, and 2.0-8.0 keV, respectively. The image 
is adaptively smoothed with a minimum signal-to-noise of 3. 



& Karachentsev 2000) late- type barred SAB(rs)bc galaxy 
with strong circumnuclear star formation. NGC 2903 
was first classified as having a central concentration of 
multiple bright hotspots superposed on a fainter back- 
ground in an early galaxy classification scheme (Morgan 
1958) based on central light morphologies (see also Sersic 
& Pastoriza 1965, 1967). These hotspots have since been 
identified as young massive star clusters in optical and 
near-IR images interspersed with numerous bright Hn 
regions (Planesas et al. 1997; Alonso-Herrero et al. 2001; 
Hagele et al. 2009) within an ^650 pc diameter circum- 
nuclear star-forming region. The star formation rate 
(SFR) in this region alone is ^0.7 M Q yr _1 (Alonso- 
Herrero et al. 2001; Leon et al. 2008). This nuclear re- 
gion is bright at wavelengths from X-rays (Tschoke et al. 
2003; Perez-Ramirez et al. 2010), UV (e.g., Popping et al. 
2010), optical and mid-IR, to CO (Sheth et al. 2002; 
Heifer et al. 2003), HCN (Leon et al. 2008), and radio 
(Wynn-Williams & Becklin 1985; Tsai et al. 2006). For 
the most part, the nuclear region displays a patchy mor- 
phology with little obvious symmetry. Perhaps the pre- 
cise center is best seen in high-resolution i7-band images 
(e.g., Perez-Ramirez et al. 2000; Alonso-Herrero et al. 
2001) where a two-arm nuclear spiral pattern is also ev- 
ident. 

The nuclear region of NGC 2903 is clearly being fueled 
by inflow of gas along its bar (Hernandez et al. 2005; Leon 
et al. 2008). The total SFR along the bar is comparable 
to the rate in the nuclear region, ^0.7 M yr _1 (Leon 
et al. 2008; Popping et al. 2010). The stellar bar extends 
roughly 1' (^2.5 kpc) from the nucleus at a P.A. of 24° 
(Sheth et al. 2002). A strong, nearly unbroken dust lane 
along the bar is clearly visible in molecular bands (Sheth 



et al. 2002; Leon et al. 2008). The dust lane leads the stel- 
lar bar (P.A. ^30°, rotation is counterclockwise (CCW) 
and the disk is inclined ~60°) but trails active star for- 
mation visible as a curved string of numerous bright Hn 
regions further downstream (Sheth et al. 2002; Popping 
et al. 2010; Knapen et al. 2002, see particularly Figure 1 
of Sheth et al.). These offsets are expected inside the 
corotation radius when gas and dust dissipates energy 
and loses angular momentum as it passes through the 
stellar bar (e.g., Martin & Friedli 1997). 

The bar is not as distinct at shorter wavelengths but 
instead appears patchy as does the outer flocculent spiral 
arm structure (Popping et al. 2010). The spiral pattern 
is clearest in Hi maps (e.g., Walter et al. 2008). Deep 
images reveal an Hi envelope extending to at least three 
times the optical diameter of the galaxy (Irwin et al. 
2009) indicating no recent interactions which tend to 
truncate the Hi disk (Higdon et al. 1998; Chung et al. 
2009). 

Hence, NGC 2903 is in many respects a textbook ex- 
ample of an isolated late-type barred spiral galaxy dis- 
playing a gas-rich star-forming nucleus and bar structure. 
This makes the galaxy an ideal target to study gas flow 
driven by the internal structures. Here we investigate the 
properties of the X-ray emitting hot gas and point-like 
sources in NGC 2903 revealed by a deep Chandra X- 
ray Observatory spectrophotometric image. We mapped 
X-ray emission in order to examine the correlation of 
the hot gas and point sources with sites of current star 
formation and of the various galactic structures in § 2. 
The detailed properties of the hot gas is analyzed in § 3, 
where we do confirm the emission asymmetry in the dif- 
fuse X-ray light first revealed by ROSAT/PSPC images 
(Tschoke et al. 2003) and visible in contemporary XMM- 
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Newton data (Perez-Ramirez et al. 2010). We search for 
an nuclear activity in § 4 and find that none of the point 
sources are coincident with the nucleus of the galaxy. 
We discount the possibility of a low-luminosity AGN in 
NGC 2903 (cf. Perez-Ramirez et al. 2010). We discuss 
that some of this emission likely traces a galactic wind 
or fountain extending above the disk of NGC 2903 origi- 
nating in the nuclear hotspot region and revealing active 
stellar feedback and its influence on the growth of nucleus 
in § 5. 

2. MAPPING THE X-RAY EMISSION FROM NGC 2903 

NGC 2903 was observed with Chandra using the ACIS- 
S instrument in imaging mode on 7 March 2010 (ObsID 
11260). The level 1 event list was reprocessed using the 
CIAO (version 4.3) tool acis-process-events to apply the 
CTI correction and to adjust for the time-dependent gain 
(with CALDB 4.3.0). Then, bad and hot pixels as well 
as bad status bits and grades were filtered out to cre- 
ate a level 2 event file. In order to identify any back- 
ground flares, a lightcurve, binned to 500 s increments, 
was created from events in the energy range of 10—12 keV 
where the telescope mirrors were insensitive to X-rays. 
Background flares were identified and removed using the 
3cr clipping method 1 which resulted in 92 ks of the final 
GTI for this observation. The angular resolution of the 
Chandra mirrors is slightly better than the pixel resolu- 
tion of the ACIS detector. Hence, we created a subpixel 
resolution image in the nuclear region by applying the 
sub-pixel event-repositioning algorithm "EDSER" of Li 
et al. (2004), which became available in CIAO v. 4. 3. For 
the rest of the galaxy, the standard 0'.'49 pixel resolution 
image was created. 

X-ray emission in star-forming galaxies like NGC 2903 
arises from various sources including bright X-ray bina- 
ries (appearing as relatively hard point-like sources), dif- 
fuse hot gas from supernovae and massive star winds, and 
underlying emission from numerous fainter unresolved X- 
ray binaries. Figure 1 displays a 4'x6' 3-band optical 
image of NGC 2903 obtained from the Sloan Digital Sky 
Survey (York et al. 2000) along with a 3-color smoothed 
X-ray image of the central 2'x3' region. The X-ray im- 
age shows numerous point-like sources as well as high 
surface brightness relatively hard diffuse emission in the 
nuclear region and more extended and softer diffuse emis- 
sion throughout much of this central region of NGC 2903. 
In this section, we first characterize the bright detected 
point sources and then describe the morphology of the 
diffuse emission and create maps of the thermodynamic 
state of the hot gas from which we identify larger zones 
that we analyze spectroscopically in § 3. 

2.1. X-Ray Point Sources 

The source finding tool in lextrct (Tcnnant 2006) was 
applied in the energy range of 0.5—8.0 keV to detect 
point sources inside the -D25 isophote. A total of 92 point- 
like sources were detected with a signal-to-noise ratio 
(S/N) above 2.4 (see Tennant 2006) and with a minimum 
of 5 counts above the background uncertainty. These 
sources are tabulated in the Appendix where we also pro- 
vide spectral analysis of the brightest sources. A spec- 
trum obtained by co-adding events from all sources was 

1 http://cxc.harvard.edu/ciao/thrcads/flarc/ 



also examined. Its rough spectral shape is that of an ab- 
sorbed power law with spectral index T = 1.95 ± 0.02 al- 
though there are some large fit residuals in the ^1—2 keV 
range likely due to contributions from (point-like) ther- 
mal emission sources. The total luminosity in point 
sources is 1.9xl0 40 erg s _1 in the 0.5—8.0 keV range. 

Visual inspection of the point source locations over- 
laid on the GALEX UV, SDSS optical, Spitzer near-IR, 
and the BIMA CO images revealed very few correlations 
with features visible at these other wavelengths. The 
notable exceptions are (1) the nuclear region which con- 
tains 8 point-like X-ray sources within a 15" (650 pc) 
radius (this region was unresolved in previous XMM- 
Newton and ROSAT X-ray observations) and (2) a 
bright object north of the nucleus, designated CX- 
OUJ093209. 7+213106, coincident with strong emission 
from the UV to mid-IR bands. 

We note that there is also strong, relatively hard, dif- 
fuse emission extending throughout the nuclear region. 
Some detected point sources in the nuclear region may be 
knots of hot gas rather than X-ray binaries traditionally 
associated with point-like sources. We inspected the X- 
ray colors and sizes of these sources to help differentiate 
compact binaries from concentrations of hot gas follow- 
ing the simple methods in Swartz et al. (2006). For the 
latter, we applied two-dimensional Gaussian models to 
the distribution of X-ray events to determine if sources 
are point-like indicating compact binaries or more ex- 
tended hot gas concentrations. None of the 8 point-like 
sources in the nuclear region were conclusively identified 
as having a thermal origin. Further analysis of the point 
sources and diffuse emission in the central 15" region is 
deferred to § 4. 

2.2. Diffuse X-Ray Emission 

We confine our analysis of the diffuse X-ray emission in 
NGC 2903 to the 2'x3' region shown in the right panel 
of Figure 1. There is no evidence of diffuse emission 
beyond this region in either the Chandra or the recent 
XMM-Newton observations. A point-source-free diffuse 
emission map was created by removing the detected point 
sources and then filling the excluded regions by sampling 
the Poisson distribution of the mean levels of their local 
backgrounds using the CIAO tool dmfilth. 

Figure 2 shows this diffuse emission map in the soft 
(0.5—2.0 keV) X-ray band. This emission map has been 
smoothed with the CIAO tool aconvolve using a Gaussian 
smoothing function with a width of 10 pixels (~5"). Also 
shown for comparison is an unsmoothed image binned 
by 4 pixels (~2") which is the pixel scale used below for 
building maps of the thermodynamic state of the gas. 
Diffuse emission is clearly present in the central region. 
The low surface brightness diffuse emission extends to a 
few minutes of arc, especially to the north of the nucleus. 

2.2.1. Galactic Structures and X-Ray Morphology 

Figure 3 shows the same 2' x 3' region as Figure 2 in dif- 
ferent wavelength images: HST/ ACS F814W, GALEX 
FUV, Spitzer 24 fim, and BIMA SONG CO. The con- 
tours of the smoothed X-ray diffuse emission from Fig- 
ure 2 are overlaid on each panel of Figure 3 to compare 
the diffuse X-ray emission to the galactic structures and 
star-forming regions visible at these other wavelengths. 
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The brightest feature at all wavelengths is the nu- 
clear region. The bar is also clearly represented in the 
HST/ACS F814W image tracing the stellar emission 
and in the CO and 24 am. images tracing cold molec- 
ular gas and dust. Popping et al. (2010) and Sheth 
et al. (2002) report an offset between these stellar and 
cold gas/dust bars and the Ha emission tracing young 
high-mass stars with the Ha emission leading the CCW- 
rotating bars. The UV image, tracing stars younger than 
~100 Myr, appears much more spotty and does not fol- 
low the bar structure. Popping et al. (2010) also note the 
UV emission is patchy and does not correlate well with 
either the Ha or 24 fim tracers of current star forma- 
tion. Instead, they find extremely young ages (<10 Myr) 
for UV-emitting clusters along the bar and much older 
(>100 Myr) for other clusters in the field. At the longer 
wavelengths, the emission extends beyond the ends of 
the bar into the trailing spiral arms and forms a distinct, 
though patchy, S-shaped pattern. 

The highest surface brightness X-ray emission also oc- 
curs in the central region of NGC 2903. There is also 
weaker emission extending beyond the nucleus in the X- 
ray image but, except for a few bright knots, the emis- 
sion is not confined to the stellar or molecular bar nor 
the S-shaped pattern visible in the IR and CO bands. 
Instead, this low surface brightness emission extends to 
the northwest of the nucleus as if filling the inner por- 
tion of the upper half of the S-shaped region of star for- 
mation traced by the 24 /mi emission. Contrary to this, 
there is no bar- like morphology along the southern half of 
the bar nor extended emission in the southwest S-shaped 
region in X-ray. We considered a difference in the col- 
umn density might account for this asymmetry; however, 
no large difference in extinction between the north and 
south halves of the bar is reported from Ha or 24 /im 
analysis (Popping et al. 2010). In summary, besides the 
nuclear region and a few isolated extended knots which 
may be associated with bright star-forming regions, there 
are no strong correlations between the X-ray morphology 
and the galactic structures seen at various other wave- 
lengths. 

2.2.2. X-Ray Thermodynamics Maps 

The lack of strong correlations between diffuse X-ray 
emission surface brightness and galactic structures visi- 
ble at other wavelengths prompted us to examine other 
properties of the diffuse X-ray emission in greater detail. 
To this end, we created maps of the thermodynamic state 
of this hot gas to examine if we see any particular mor- 
phological features in the X-ray properties. 

Specifically, we created a point-source-free image 
binned by 4 pixels (~2"). Then, we extracted a spectrum 
from a circular region around each pixel in the binned 
image. This procedure follows the method described 
in Randall et al. (2008), who employed the technique 
of O'Sullivan et al. (2005) and Maughan et al. (2006). 
The radius of each region was chosen to contain at least 
500 counts after background subtraction. A region to 
the northeast outside of the D25 annulus was used for 
the background. We constrained the maximum radius of 
each region to be <30"; we discard regions needing >30" 
to contain >500 net counts. Overall, the resulting range 
of radii is 5"— 30". This is much larger than the pixel 
scale representing the resulting parameter values; this 



is in a sense equivalent to "smoothing" the parameter 
maps. 

Each spectrum is then fitted using the energy range 
between 0.5—2.0 keV, because there are very few counts 
detected above 2 keV after the background subtraction. 
We fitted with an absorbed collisionally ionized optically 
thin plasma model, namely phabs* mekal in XSPEC. 2 
Due to low counts in each spectrum, we use Poisson 
statistics to fit spectra and we fixed the metallicity to 
the solar value following the XMM-Newton Reflection 
Grating Spectrometer (RGS) analysis of Perez-Ramirez 
et al. (2010); see also Pilyugin et al. (2004) and refer- 
ences therein. Applying the same model to the spectra 
of each spatial region allows us to spatially resolve rela- 
tive changes in the thermodynamic properties of the gas. 

Figure 4 shows the resulting X-ray temperature, kT e , 
absorbing column density, Ah, electron density, n e , and 
pressure, P, maps of the hot gas over the same regions 
as displayed in Figure 2. We also overlay the HST/ ACS 
F814W contours in order to visualize the underlying 
galactic structures. 

The measured temperature shown in Figure 4 ranges 
from 0.2 to 0.3 keV with errors of 0.03 keV on average. 
This is a typical temperature for hot gas (0.1—0.7 keV) in 
star- forming regions in normal or starburst galaxies (e.g., 
Tyler et al. 2004; Strickland et al. 2004). The nuclear 
starburst region has relatively higher kT e than the sur- 
rounding areas. The highest temperature in the nuclear 
region (yellow in the figure) is elongated in a roughly 
NW to SE direction nearly perpendicular to the bar and 
extending to the NW. This extension to the northwest is 
also seen in the XMM-Newton image presented in Perez- 
Ramirez et al. (2010). However, there is no large differ- 
ence in temperature among these regions overall, which 
is also consistent with what has previously been found 
in the other star- forming galaxies (e.g., Swartz et al. 
2006; Wang et al. 2009). Interestingly, the temperature 
along the (northern) bar decreases whereas the surface 
brightness remains roughly constant. The temperature 
is higher to the northwest of the northern bar (within 
the S-shape). 

The measured Ah is in the range of (2—8) x 10 21 cm~ 2 , 
which is much larger than the Galactic Ah value of 
3xl0 20 cm' 2 (Dickey et al. 1990) along the line of sight. 
In fact, this is larger than the intrinsic Ah value for 
NGC 2903 of -5 x 10 20 cm" 2 derived from the Hi Nearby 
Galaxy Survey data (Walter et al. 2008) estimated us- 
ing their Equation (5). We note that applying a single 
temperature model to the multi-temperature gas seen in 
galaxies often tends to result in a higher column density 
because of the need to suppress low-energy emission in 
the model (e.g., Yukita et al. 2010). Due to the low fitted 
model temperatures, the Ah uncertainties are somewhat 
large, about ±(0.5— 1.0) xlO 21 cm' 2 . The nuclear region 
and the dust lane to the east of the nucleus have the 
highest Ah- The north side of the bar also shows a rela- 
tively higher column density, and the regions to the west 
of the nucleus, where fewer and weaker dust lanes are vis- 
ible in the optical Hubble Space Telescope (HST) image, 

2 The mekal model was used here because the solution is com- 
puted "on-the-fly" as opposed to being pre-defined on a finite model 
parameter grid as does the apec model so that it does not suffer 
from a discretization problem. 
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Fig. 2. — Left: the central 2'x3' of the point source-free X-ray soft (0.5—2.0 keV) diffuse emission map. The image was smoothed 
using the CIAO tool, aconvolve with Gaussian a width of 5". The contours are also plotted in square root scale, specifically 0.075, 0.132, 
0.307,0.59, 1.0 and 1.52 count pixel. The pixel scale of this image is 0''492. Right: the same as the left image, besides binned by 2" scale 
rather than smoothed. The color bars for the both images indicate the square root of the intensity with the unit of count. 



have relatively lower A^h values. However, the spectral 
properties, both Ah and kT, along the south side of bar 
does not show a correlation with the underlying galactic 
structures, similar to the lack of correlation with these 
features in the X-ray surface brightness. 

The X-ray electron density, n e , is estimated from the 
model fit parameters in the following way. The spectral 
model normalization, K, is related to the emission inte- 
gral, J n e nndV ~ n^fV, where V is the X-ray emitting 
volume and / is a volume filling factor. This volume 
filling factor represents the fraction of the volume actu- 
ally occupied by X-ray emitting gas. We assume that 
the number densities of ions, rii, and electrons are equal. 
This is a reasonable assumption for hydrogen dominated 
plasma. The electron number density, then, can be esti- 
mated as rig oc K/(fV). 

In order to derive an electron density map, we assume 
that the hot gas fully occupies (/ = 1) a cylindrical 
volume represented by a scale height, h = 200 pc, of 
the disk. The electron density then scales linearly with 
This is, of course, not the only possible scenario; 
the gas (or some fraction of the gas) may extend above 
the disk in the form of an outflow. This possibility is 
examined further in § 5. 

Using the derived n e and the flux-weighted mean X- 
ray temperature T e , the hot gas pressure P/k — 2n e T e 
can be derived. The pressure is also highest in the nu- 
clear region, approaching 25 xlO 5 K cm~ 3 , and is more 
centrally-peaked than either the temperature or the elec- 
tron density alone. The pressure distribution is smoother 
beyond the nuclear region approaching a roughly con- 
stant value of (5—7) xlO 5 K cm~ 3 . 

In summary, Figure 4 shows that the highest density, 
temperature, and pressure region is the nuclear starburst 
region. This is also the highest surface brightness region 



(Figure 2). The density is relatively high in the region 
along and to the northwest of the northern portion of 
the stellar bar. We do not see an analogous feature cor- 
responding to the south side of the bar. Lastly, the pres- 
sure map indicates that only the nuclear starburst region 
is highly pressurized if we assume an equal scale height 
for all the regions. In fact, it is possible that the hot 
gas is flowing out into a low density halo region, which 
would imply a different morphology of the gas volume 
than considered here. 

We point out that the fitted column density and tem- 
perature parameters are often anti-correlated. There are 
low column density (blue) regions, where the temper- 
ature is high (yellow and orange regions west of the 
nucleus). Also, the temperature is lower (red) where 
the column density is higher (yellow). Another, though 
weaker, positive correlation sometimes occurs between 
column density and electron density (proportional to the 
square root of the model normalization) in these single- 
temperature models. We show the column density of all 
the pixels in the map against their corresponding tem- 
perature to examine this behavior in Figure 5. 

The left panel of Figure 5 clearly shows that most re- 
gions follow the anti-correlation with higher Ah values 
tending to correspond to lower temperatures. These re- 
gions are located in the bottom left side of the figure. 
There are also regions that do not follow this trend, 
spreading into the top right of the plot. We selected 
three groups in this parameter space as indicated by 
lines in the left panel of Figure 5 based on the group- 
ing of the regions apparent to the eye. The right panel 
of Figure 5 shows the spatial map of the corresponding 
groups defined in this way. As anticipated, the nuclear 
region corresponds to the top right group consisting of 
the hottest, most heavily absorbed gas with clearly dif- 
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Fig. 3.— Central 2'x3' of the NGC 2903 of FST/ACS F814W (top left), GALEX FUV (top right), CO (bottom right), and Spitzer 24 /im 
(bottom right) images. The color bars indicate a square root of intensities. The HST image has the pixel scale of 0'.'05 and the intensity 
unit of electron; le corresponds to a fluence of 7.1xl0 — 20 erg cm -2 A -1 , that is a specific flux F\= l.OxlO" 22 erg s — 1 cm -2 A" 1 . The 
pixel scale of the GALEX FUV image is 1'.'5 with the intensity units of counts s — 1 . The BIMA SONG CO image has the pixel scale of 1". 
The unit of image is Jy (beam) -1 km s — 1 . The Spitzer 24 fim image has the pixel scale of l'/5, and the map intensity is in Jy sr _1 . The 
contours depict the intensity of the smoothed X-ray diffuse emission shown the left panel of Figure 2. 
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Fig. 4. — Same region as in Figure 2 of the temperature (top left), hydrogen column density (top right), electron density (bottom left), 
and pressure (bottom right) of the hot gas derived by applying an absorbed single temperature model (see the text). The units for each map 
are in keV, 10 22 cm -2 , cm -3 , and 10 s K cm -3 , respectively, as denoted by the color bar at the top of each panel. One pixel corresponds 
to ~2", which is the same as the right panel of Figure 2. However, the spectral extraction regions used to define the pixel values are larger 
than 2", see the text for details. The regions which did not contain enough counts for spectral fitting appear black. We note that a single 
temperature model can result in higher TVjj and normalization, therefore resulting in a higher density. In order to compare with the galactic 
structures, the HST/ ACS F814W contours are overlaid in green. 
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Fig. 5. — Left: The iVjj vs. kT e plot for the fits of an absorbed single temperature model using 0.5—2.0 keV. In general, there is an 
anti-correlation between jV H and kT e , as a higher JVjj results in a lower kT e value. We divided the parameter space into 3 groups as 
indicated with straight lines. Right: the three parameter groups defined in the left panel arc mapped in spatial coordinates. The light gray 
zone indicates the right top group, and the dark gray zone depicts the bottom left group, and the black zone is the middle group in the 
left panel. These two plots show that the nuclear region has both high Nn and kT e values that do not follow the correlation seen in the 
rest of region. 



ferent properties than the other two regions. The middle 
and bottom left groups also show strong spatial correla- 
tions. The middle group regions are located to the east 
and the bottom left group are located primarily to the 
west of the nuclear region. 

Although these regions were found using solar abun- 
dances, we note that similar structures are evident when 
changing abundances to a different value, such as Z = 
0.05 Zq. In order to understand why these regions are 
different, we analyze the X-ray spectra extracted from 
these three groupings, or zones, where we are able to 
apply more complex models in the next section. 

3. FURTHER SPECTRAL ANALYSIS OF THE DIFFUSE 
EMISSION 

We have shown in the previous section, using detailed 
maps of the thermodynamic state of the hot, X-ray- 
emitting gas, that the spatial complexity observed at 
other wavelengths is clearly present at X-ray energies. 
On the other hand, there are no strong correlations be- 
tween the X-ray morphology and the underlying galac- 
tic structure with the notable exception of the nucleus. 
There is a general trend of lower temperature and lower 
column density to the west of the nuclear region relative 
to east of the nuclear region. In this section, we will 
analyze in greater detail the spectral properties of three 
zones: the ~15" nuclear starburst region, and the two 
extended diffuse emission regions to the east and to the 
west of the nucleus. 

3.1. The Nuclear Zone 

We performed X-ray spectral analysis of the diffuse nu- 
clear emission defined as the light gray zone of Figure 5 
from the JVh — kT e correlation. This zone is roughly 15" 
(650 pc) in radius. We first fit the same absorbed ther- 
mal model using the same energy range (0.5—2.0 keV) 



as before to validate the results in our thermodynamic 
maps. We obtained Nn =5.8lg'3Xl0 21 cm~ 2 , and 

fcT e =0.22to;oi keV, which is consistent with the results 
from the maps. However, with many more counts in the 
source spectrum, this simple absorbed thermal model 
with a solar abundance is not well fitted to the data 
(Xl ~ 4). 

In fact, there are enough counts to extend the model 
fit to higher energies, 0.5—6.0 keV, resulting in 2205 ±48 
net counts after background subtraction. The same local 
background, taken from outside the -D25 isophote, that 
was used for creating the thermodynamic map was used 
here. The spectrum was grouped to contain at least 40 
counts in each spectral bin, and the \ 2 statistic was used 
for determining the goodness of fit. 

Including the higher energy photons resulted in raising 
the temperature (kT e =0.5) and lowering the column den- 
sity (2xl0 21 cm~ 2 ) of the best-fitting single-temperature 
model but the fit is unacceptable (x 2 ~ 6). Next, we 
allowed the elemental abundances in the mekal model 
component to vary keeping their ratios fixed to the solar 
ratios given by Anders & Grevesse (1989). This dramati- 
cally improved the fit to xl ~l-9 with one fewer degree of 
freedom and resulted in raising the best-fitting tempera- 
ture to fcT e =0.58 keV and lowering iV H to 3xl0 20 cm -2 . 

The resulting abundances were subsolar, Z = 0.05 Zq, 
which is inconsistent with the XMM-Newton/RGS re- 
sult (Perez-Ramirez et al. 2010) and is problematic if 
the hot gas originates from a-element-rich core-collapse 
supernovae as expected. On the other hand, the widely- 
applied Weaver model (Castor et al. 1975; Weaver et al. 
1977) of hot bubble formation predicts the majority of 
hot gas is mass-loaded from the surrounding cold clouds 
rather than coming from the metal-rich shocked super- 
novae gas. It is also known that a simplistic model some- 
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times gives very low abundances when applied to the 
relatively low CCD resolution imaging spectra (Weaver 
et al. 2000). Hence, we do not argue strongly that the 
resultant abundance estimate is realistic and decided to 
fix abundances at the solar value, consistent with Perez- 
Ramirez et al. (2010), in the remainder of our analysis. 

The fitting residuals indicate that the single- 
temperature model best fits the data only between 0.7 
and 1.2 keV. Therefore, we added a second model compo- 
nent, either a power-law component or a second thermal 
model, to account for the softer and harder emission. 

Adding a power-law model resulted in a significantly 
improved fit statistic (xt The power-law index 

is r=2.65. This value is similar to that found in the 
nuclear starburst regions of NGC 1365 (Wang et al. 
2009) and NGC 5135 (Levenson et al. 2004). These au- 
thors interpreted this hard emission as unresolved point 
sources; however, the power-law slope is much steeper 
than the slope of the detected point sources in NGC 2903 
(r ^1.95; see the Appendix). Also, fitting only the data 
in the range 2 and 6 keV gives an even steeper index of 
~4. We argue that the harder X-ray emission is unlikely 
to come from unresolved point sources and is instead a 
hotter thermal component. 

Adding a second thermal component instead of a power 
law results in a xt ~2, but not yet an acceptable fit. 
The two temperature model results in a second temper- 
ature of ~3 keV, which may be interpreted as shocked 
gas directly associated with a "wind fluid" as observed 
in powerful starbursts like M82 (Strickland & Heckman 
2007, 2009). This is a more realistic interpretation than 
the unresolved point source (power law) interpretation 
of the hard emission for the nuclear region of NGC 2903. 
To obtain a better fit statistic, however, a third thermal 
component is required. The resultant fitting parameters 
are tabulated in Table 1 where the errors are quoted at 
the la (68% confidence level). Figure 6 shows the spec- 
trum of the nuclear region for this 3-temperature model 
illustrating the complex line structures due to O and Fe 
emission between 0.6—1.2 keV, Mg at 1.3 keV and Si at 
1.8 keV as well as a hard continuum above 2 keV. Note 
the fit residuals indicate some of these line features are 
still poorly reproduced by this model with abundances 
fixed to their solar values. The two lower temperatures 
for this model are ~0.2 and 0.6 keV, which are typical 
values for star-forming regions and the halos of starburst 
galaxies (Strickland et al. 2004; Yukita et al. 2010). This 
model results in a column density Ah=1.4x10 21 cm~ 2 . 

We also derived other physical properties of the hot 
gas from the spectral fits; namely, the electron density 
n e and pressure P, derived in the same manner as in 
§ 2.2.2. Again, we note that the assumed geometry and 
hence hot gas volume can introduce large uncertainties 
on these derived properties. In addition to these two 
properties, we derived the hot gas mass, Mx —n e ^,m p V, 
where m p is the proton mass and /i=1.4 is the mass per 
proton, the thermal energy £7 t h = 3n e VkT e , and the 
cooling time r=E t h/ '-^boi- Here, Lbol is estimated from 
the absorption-corrected X-ray hot gas luminosity over 
the energy range between 0.01— 20 keV. The results are 
listed in Table 2. It has been assumed that the individ- 
ual temperature components are in pressure equilibrium 
in order to estimate the electron density and density- 
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Fig. 6. — Spectrum of the nuclear region defined by ./Vh — kT e 
correlation of Figure 5. The spectrum was grouped to achieved 
at least 40 counts in each bin. A background spectrum, extracted 
from a region to the northeast of the galaxy and outside of the D25 
isophote, has been subtracted. An absorbed three temperature 
component mekal model with solar metal abundances is shown as 
the solid curve. The lower panel shows the contributions to the \ 2 
fit statistic. 

dependent properties. 

Finally, the X-ray spectral properties of the central re- 
gion of NGC 2903 can be compared to the properties 
previously reported by Perez-Ramirez et al. (2010) us- 
ing a deep XMM-Newton observation. XMM-Newton 
cannot resolve the point sources. Therefore, we fit- 
ted the Chandra spectrum of the central 15" radius 
region including point sources. The best-fit model to 
the XMM-Newton data is of an absorbed power law 
plus thermal component, phabs (powerlaw+mekal), with 
power law T=2.1, kT e = 0.5 keV, and I>x(o.5-8 keV) — 
2.3xl0 39 erg s _1 , (Perez-Ramirez et al. 2010). The 
Chandra spectrum has higher luminosity, £x(o.5-8 keV) 
= 3.2xl0 39 erg s _1 and is characterized better with a 
combination of absorbed disk blackbody and thermal 
models, phabs (diskbb+mekal) , with kT- m =1.3 keV and 
kT e = 0.5 keV. When a power-law model is used instead of 
a disk blackbody, the fitted power law index is lower than 
for the XMM-Newton observation, T =1.8, again driven 
by the hard X-ray spectral shape. Both spectra suggest 
the absorption column is very low, at the Galactic value, 
and that the hot gas temperature is higher than that of 
the map we derived in the previous section. Note that 
the diffuse emission only accounts for about 1/3 of the 
total luminosity in the Chandra data (cf. Table 1) with 
the remainder contributed by point-like sources. The 
largest contribution, £x(o 5-8 keV) = 1.4xl0 39 erg s _1 , 
is from CXOJ023210. 1+213008, which is also fitted well 
with a disk blackbody model (see the Appendix) with 
fcTj n =1.03 keV. Since X-ray binaries tend to dominate 
the luminous X-ray point-source population and they 
are often time variable, the difference in luminosity es- 
timated here compared to the XMM-Newton results of 
Perez-Ramirez et al. (2010) may be due simply to vari- 
ability in one or more of these bright X-ray point sources. 

3.2. Extended Diffuse X-ray Emission 

Following the results of § 2.2.2, we divided the low 
surface brightness diffuse emission surrounding the nu- 
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TABLE 1 

Fit Parameters of X-Ray Diffuse Emission 



Parameter 


Nucleus 


West 


East 


Net count 


2448±49 


1928±53 


2512±61 


Fitted energy range (keV) 


0.5-6.0 


0.5-4.0 


0.5-4.0 


n H /10 22 (cm- 2 ) 


14+ 004 


o.2otr 6 6 


n qy+0.05 
U -^'-0.05 


kT el (keV) 


n 10+O.O2 
U - 18 -0.02 


o.is±g:8i 


1Q+ 001 


Ki/icr 4 


0.38t°,; 23 


1 50+0.10 

1 - oo -0.07 


9 =7+1-32 
z -°'-0.89 


kT e2 (keV) 


u ' oo -0.02 


u - oo -0.04 


56+ (U)3 


K2/IO- 4 


0.36I 007 


0.28t°;° 5 5 


u -**-0.06 


kT e3 (keV) 










S0+ 03 






X 2 /dof 


60.1/42 


59.2/41 


75.36/49 


L X0 5-8 o/ 1038 ( Cr S O 


10.8+ij 


7 -+0.0 
'•'-1.4 


9.81°;° 


^kt.s-s.o/io 38 ( cr e O 


17.9 


140 


194 



TABLE 2 

Properties of X-Ray Diffuse Emission 



Parameter Nucleus West East 



Volume [/] (10 63 cm 3 ) 9.6 56 67 

M/~ 1/2 1 (cm- 3 ) 1.07, 0.36,0.06 0.08,0.03 0.10,0.03 

P/k[f' 1/2 ] (10 s K cm- 3 ) 46.9 3.5 4.3 

Eth[f +1/2 ] (10 53 erg) 93.5 40.4 57.2 

M x [f +1/2 ] (10 5 M Q ) 7.6 32.1 44.5 

r (Myr) 63.2 6.6 13.0 



clear zone into two regions, namely the (roughly) east 
and west zones based on the ./Vh — kT e correlation map 
shown in Figure 5. These zones do not show strong cor- 
relations with the underlying galactic structures seen at 
other wavelengths. However, the locations of dust lanes 
relative to the (optical) starlight from the spiral arms 
(e.g., the SDSS composite image of Figure 1) suggest 
the eastern region lies closer along our line of sight than 
the western region. Importantly, if some of the extended 
(non-nuclear) hot gas diffuse emission lies above (below) 
the plane of the galaxy, then we would expect the rela- 
tive absorbing column density to be lower (higher) in the 
west (east) zone. This is indeed what is deduced from 



the thermodynamic maps used to define these two zones 
(Figure 5). 

In order to compare to the hot gas in the nuclear star- 
burst region, the same sequence of models were applied 
to the hot gas spectra in the east and west zones as was 
done for the nuclear region with abundances of the ther- 
mal component fixed to solar. 

We first applied the same single-temperature absorbed 
thermal model as used for the thermodynamic map over 
the same energy range. The spectra were grouped to have 
at least 40 counts in each bin, and x 2 statistic is used. 
The same background spectrum that was used above was 
applied. 

Again, we have reproduced consistent results, specif- 
ically A H = (4.6±0.03) x 10 21 cm- 2 and kT e = 
0.21±0.01 keV for the east region, and Ah = 
(4.3±0.03)xl0 21 cm" 2 and kT e = 0.19±0.01 keV for the 
west region. Likewise, the fitting statistics are not ac- 
ceptable (xf, >2) with larger number counts in the source 
spectra. Again, there is some emission above 2 keV, 
hence we expanded the energy range to 0.5—4.0 keV for 
subsequent spectral analysis. 

In both zones, including the higher energy range data 
did not change the fitting results significantly when fit- 
ting to a single-temperature thermal model with solar 
abundance. 

Adding a second component (either power law or sec- 
ond thermal component) improved the fit statistic (x 2 
~1.4) in both the zones. However, the best-fitting power 
law indices in both zones are r ~ 3.7. Such steep slopes 
are unphysical for unresolved point sources. Therefore, 
we applied two-component thermal models. The fitting 
results are listed in Tables 1 and 2 and the observed data, 
spectral model, and fit residuals are shown in Figure 7. 
Adding a third thermal model component did not im- 
prove the fit according to the F-test. 

The best-fit temperatures are ~0.2 keV and ~0.6 keV, 
which are consistent with the two lower temperatures 
in the nuclear zone and with typical star-forming re- 
gions and the halos of starburst galaxies (Strickland et al. 
2004; Yukita et al. 2010). The dominant contribution is 
from the lower-temperature component, which is consis- 
tent with the single-temperature model. Adding a sec- 
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ond temperature component lowers the absorbing col- 
umn densities in the two zones with in the eastern 
zone slightly higher than in the western zone though they 
agree within errors. 

4. THE NUCLEAR STARBURST AND A SEARCH FOR AN 
ACTIVE NUCLEUS 

Images of the nuclear region in several wavelengths are 
shown in Figure 8. Each image is 30"x30" centered 
at R.A.=9 /l 32 m 10!l, dccl =21°30'03". The emission is 
rather patchy at most wavelengths but a rough ring-like 
structure, ~650 pc in diameter, has been reported in the 
literature (e.g., Alonso-Herrero et al. 2001). Although 
the size of the ring is debated as one-half this size even 
after deprojection (Comeron et al. 2010). The 24 fim 
image has the lowest resolution, ~5". The 24 /xm flux 
is dominated by the star-forming regions to the north of 
the center. The 8 /im image shows more of a ring-like 
morphology and is similar to comparable-resolution radio 
images (not shown; see, e.g., Wynn- Williams & Bcck- 
lin 1985). The JZST/NICOMS 1.6 and HST/ ACS 
F814W images show the numerous star clusters present. 
The F814W image also reveals that there are complex 
dust lanes in the east side of the nuclear starburst re- 
gion. The Ha image also reveals the ring-like morphol- 
ogy (Alonso-Herrero et al. 2001) consisting of many Hn 
regions. 

The high resolution Chandra image of the nuclear re- 
gion reveals that there are a number of point-like sources 
embedded in the soft X-ray diffuse emission. Most do 
not correlate spatially with features visible at compara- 
ble resolution at other wavelengths as expected if these 
sources are X-ray binaries. Overall, the X-ray morphol- 
ogy does not trace the ring-like shape but is center-filled 
and circular rather than the elliptical shapes seen at 
other wavelengths that roughly match the inclination of 
the galaxy. The soft X-ray diffuse emission seems to be 
affected by the dust lanes, showing less emission from the 
east side of the center and more absorption, appearing 
green in the tri-color image on the east edge of the star- 
burst region where the strong dust lane can been seen. 

As Figure 8 shows, the mass center of NGC 2903 is 
not well defined because of complex structures and the 
lack of an obvious active galactic signature. NASA/IPAC 
Extragalactic Database lists the galactic center as 
R.A.=9 /l 32 m 10!ll, decl.=21°30'03'/0, which is obtained 
from the Two Micron All Sky Survey. Trachternach 
et al. (2008) estimated the galactic center position by fit- 
ting tiltcd-ring models to velocity fields derived from Hi 
observations (R.A.=9' l 32 m 10!0, decl.=21°30'0275), by 
fitting ellipses to the surface brightness of the Spitzer 
3.6 urn image (R.A.=9 ,l 32 m lffil, decl.=21°30'04'. / 9), and 
by the location of a point source detected in the radio 
continuum (R.A.=9' l 32 m 10?l, decl.=21°30'04'/3). We 
carefully registered the HST/ ACS image to the 2MASS 
coordinates but found none of the bright star clusters or 
ionized regions arc coincident with any of these four po- 
sitions. No individual X-ray point sources are detected 
at any of the four positions either. 

We place an upper limit X-ray luminosity for any un- 
detected source in the nuclear region using 1" source 
extraction radii. The extinction toward the nuclear 
region is relatively moderate, reportedly ranging from 
Ay ^2 — 4 mag (Alonso-Herrero et al. 2001; Quillen 



& Yukita 2001) for the Hn regions to Ay ~ 5.8 mag 
(Alonso-Herrero et al. 2001) for the average over the 
full 14" diameter nuclear region corresponding to hy- 
drogen column densities of (4.6— 11.2) x 10 21 cm -2 based 
on the standard Galactic gas-to-dust ratio Nn/Ay = 
1.9xl0 21 cm" 2 mag^ 1 (Maiolino et al. 2001, and ref- 
erences therein). (Note these values are considerably 
higher than our best-fit values derived from fits to the 
X-ray spectra of the nuclear region.) A recent high reso- 
lution CO observation indicates that the column density 
can be as high as (6±3)xl0 22 cm -2 (P.-Y. Hsieh, private 
communication, 2011). Therefore, we assume a column 
density of 10 22 cm~ 2 to estimate a conservative upper 
limit X-ray luminosity from an undetected point source. 
Assuming a power law spectral index of T=1.7, the 90% 
confidence upper limit luminosity from our Chandra ob- 
servation is ^1.5xl0 38 erg s _1 . This low value is consis- 
tent with the lack of any evidence at other wavelengths 
for an active nucleus. 

5. DISCUSSION 

The X-ray emission observed by Chandra from within 
the central few minutes of arc of NGC 2903 is com- 
posed of both hard-spectrum point-like sources (observed 
L x ~ 1.9 x 10 40 erg s" 1 in the 0.5-8.0 keV range) and 
soft and hard diffuse extended emission (total observed 
Lx ~ 0.3 x 10 40 erg s" 1 ). The intrinsic luminosity of 
the detected point sources is about 25% larger based on 
our spectral fits but the intrinsic luminosity of the soft 
diffuse component can be a factor ^20 higher (Table 1) 
because the observed spectrum is so soft. The observed 
diffuse X-ray luminosity from the central 650 pc region 
alone is ^10 39 erg s _1 with about 30% of this amount 
contributed by a hot component with kT e ~ 3.6 keV. 

Much of the soft diffuse component is likely created 
by the collective effects of massive star winds and su- 
pernovae in star-forming regions located throughout the 
disk of NGC 2903. We discuss below whether or not the 
hard diffuse component may be due to a concentrated 
hot wind fluid emanating from the central starburst re- 
gion, the possibility of outflow from this region into the 
halo, and whether or not some portion of the soft diffuse 
component originates from this outflow. 

5.1. Hard Diffuse X-Ray Emission 

Strickland & Heckman (2007, 2009) have done a thor- 
ough analysis of the X-ray emission from the central 
region in the nearby starburst galaxy, M82. We find 
many parallels between the X-ray properties of the cen- 
tral regions of NGC 2903 and those of M82, as given by 
Strickland & Heckman, with the notable exception that 
NGC 2903 does not so readily display copious galactic 
winds as does M82. Part of the reason is that NGC 2903 
is viewed at an inclination of 60° compared to the edge-on 
orientation of M82 but the wind power (or, equivalently, 
mass loss rate) is apparently much weaker in NGC 2903 
compared to M82. 

The hard X-ray diffuse emission in M82 is confined to 
the nuclear starburst region (r ~ 500 pc), whereas the 
soft X-ray diffuse emission extends into the extraplanar 
regions above and below the stellar disk. A similar dis- 
tribution can be seen in NGC 2903. The hard X-ray 
diffuse emission is only seen in the central ^$15" (650 pc) 
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Fig. 8. — Circumnuclear region (30" X 30") of NGC 2903. North is up, and east is to the left. Top: (from left to right) the Spitzcr 
24 /im, Spitzer 8 fim, FST/NICMOS 1.6 /jm HST/ ACS F814W. Bottom: HST/ ACS continuum-subtracted Ha, XMM-Newton OM UV, 
Chandra 0.5—8.0 keV X-ray subpixel image, and energy color-coded Chandra X-ray image. The X-ray colors indicate the 0.5—1.0 keV 
(red), 1.0—2.0 keV (green) and 2.0—8.0 keV (blue) energy bands. Red circles indicate the positions of the detected point-like X-ray sources. 
The green circle depicts the position of the galaxy center as listed in the NASA/IPAC Extragalactic Database. The radius of a circle is 
l'/5. 



radius nuclear region while the soft X-ray emission ex- 
tends ~5 kpc or more mostly to the north and west of 
the nucleus. Strickland & Heckman (2007, 2009) also 
argue that the hard continuum emission is of thermal 
origin; however, as with NGC 2903, the hard contin- 
uum from M82 is formally better fitted with a power 
law model, so that non-thermal processes such as in- 
verse Compton scattering may contribute to the hard 
X-ray emission. There are too few net X-ray source 
counts detected from NGC 2903 to quantify any underly- 
ing non-thermal emission component or to use emission 
lines from highly-ionized metals to better constrain the 
thermal component of the hard emission as Strickland & 
Heckman (2007) were able to do for M82. 

The hard continuum emission temperature from M82 
measured by Strickland & Heckman (2007) is ^3.6 keV. 
This is the same hard-component temperature we found 
for the nuclear starburst region in NGC 2903. However, 
the central hot gas pressure deduced from our analysis of 
the X-ray spectrum is 5xl0 6 K cm -3 whereas Strickland 
& Heckman (2009) estimate pressures closer to 10 7 to 
3xl0 7 K cm~ 3 from models that reproduce the observed 
hard X-ray luminosity from M82. 

Following the notation of Strickland & Heckman 
(2009), we define the thermalization efficiency, e, as the 
fraction of mechanical energy from star formation that 
is imparted to (heats) the surrounding gas (the remain- 
der is lost radiatively) and the mass-loading factor, (3, 
as the ratio of the total mass of heated gas to the mass 
directly ejected by supernovae and stellar winds. The 
central temperature in a starburst is then (see the re- 
view by Veilleux et al. 2005) 



T c = 0.4 



kpM w ' 



(1) 



where /i = 1.4 is the mean molecular weight, mu the 
proton mass, Ly/ the rate of energy input, My/ the 
rate of mass ejected from stars and supernovae, and k 
is Boltzmann's constant. Furthermore, by equating the 
rate of energy deposited in the surroundings, eLy/, to 
the asymptotic kinetic energy loss rate, (l/2)f3My/v'^ D , 
it follows that the starburst drives a wind with terminal 
velocity 

2 _ 2eL w _ 5kT c 
V °° ~ f3M w ~ MTO H ' 

Since both Ly/ and My/ scale linearly with the SFR, 
T c and Voo are independent of SFR but do depend on 
the environment through e and /3. Clearly, both quan- 
tities decrease as more cold gas mass is loaded into the 
wind fluid and as radiative losses increase. Values for 
Ly/ and My/ as a function of time can be taken di- 
rectly from the Starburst99 (Leitherer et al. 1999) stel- 
lar population synthesis models. Using our assumption 
of solar metallicity and a Salpeter initial mass function 
(Salpeter 1955) over the stellar mass range 1 — 100 M , 
and assuming a continuous rate of star formation begin- 
ning at some time t , Ly/jMy/ is approximately con- 
stant after an initial transient phase lasting <20 Myr 
so that T c ~ 2.5 x 10 8 (e//?) K. Matching this to the 
observed temperature of the hot component in the nu- 
clear starburst region of NGC 2903, (e//J) ~ 0.15 so that 
Woo ~ 3000(e//3) 1 / 2 ~ 1100 km s" 1 . For reasonable val- 
ues of e on the range 0.5-1.0 (Veilleux et al. 2005, and 
references therein) appropriate to the moderate-density 
central region of NGC 2903, the mass-loading factor is 
relatively high, 3.3 < j3 < 6.7, compared to typical values 
derived by Strickland & Heckman (2009) from numerical 
models for M82 (j3 ~ 2.5). 
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The hard X-ray luminosity scales as (f3Mw) 3 /cR+L\y 
(see Strickland & Heckman 2009, for a derivation) where 
Ri, is the radius within which the starburst energy and 
mass are injected. Taking the mean value of f3 — 5 for 
NGC 2903 and /3 = 2.5 for M82 (Strickland & Heckman 
2009), assuming e is the same for both galaxies (they 
likely differ by less than a factor of 2), and noting that 

M^/Lw oc SFR 2 , we can solve for R* for NGC 2903 
given the central SFRs and hard X-ray luminosities for 
the two galaxies and R+ = 300 pc for M82 (Strickland & 
Heckman 2009). This gives i?* ~ 780 pc for the central 
starburst in NGC 2903 which is close to our estimated 
size of the star-forming region of 650 pc. Therefore, we 
conclude that the detected hard X-ray diffuse emission 
from NGC 2903 is similar to that from M82 suggesting 
a possible outflow originates from the central starburst 
provided the current star formation activity has persisted 
for ^20 Myr or longer. This implies the hot wind fluid 
present in the nuclear region of NGC 2903 dominates 
the energetics of the flow as is implicit in the original 
Chevalier & Clegg (1985) model on which the work of 
Strickland & Heckman (2009) is based. It is natural to 
ask whether or not the hot gas carries enough energy to 
escape from the galaxy potential. The escape velocity, 
v e , is roughly 2.6—3.3 times the local circular velocity 
(e.g., Veilleux et al. 2005). The circular velocity of the 
inner 7 kpc (-6') in NGC 2903 reaches 230 km s" 1 (de 
Blok et al. 2008), yielding v e = 600-760 km s" 1 . There- 
fore, the high-temperature component of the wind fluid 
is likely to flow outward and eventually to escape from 
the galaxy disk. 

5.2. Soft Diffuse X-Ray Emission 

Unlike M82, direct detection of any galactic wind em- 
anating from the central regions of NGC 2903 is com- 
promised due to the inclination of the galaxy disk to 
our line of sight. There is soft diffuse X-ray emission 
with characteristic temperature 0.2—0.6 keV extending 
beyond the central starburst with a morphology that 
does not closely match the underlying structures seen 
at other wavelengths (§2.2.1). Much of this gas probably 
originated from supernovae and stellar winds from mas- 
sive young stars distributed throughout the stellar disk 
of NGC 2903. Its relatively soft spectrum is likely due 
to dilution of hotter gas by thermal evaporation of sur- 
rounding cold ISM as in the standard model for the for- 
mation of interstellar bubbles (Castor et al. 1975; Weaver 
et al. 1977). However, a portion of this soft X-ray emis- 
sion may come from a galactic wind, also diluted through 
mass-loading and further cooled by adiabatic expansion, 
but lying above the disk of the galaxy. 

Based on the thermodynamic maps analyzed in §2.2.2, 
the soft diffuse X-ray emission can be divided into three 
distinctive zones, the nuclear, west and east zones. We 
note that the hydrogen column density is lower in the 
west than in the east zone and lowest in the nuclear zone. 
This may suggest the geometry of the hot gas: in the nu- 
clear zone, the gas extends into the lower density halo. 
The gas in the east zone is mostly confined to within the 
disk. Because the galaxy is tilted so that the east side is 
closest to us, any emission from a galactic wind along a 
line of sight through the east zone would be attenuated 
by the disk giving the higher Nn value. Conversely, the 



TABLE 3 
Properties of X-Ray Point Sources 





R.A. 
(J2000) 


Decl. 
(J2000) 


Counts 


S/N 


LJ/10 38 
(erg s _1 ) 


1 


09 


31 


54.5 


21 


27 44.6 


22.5 


3.8 


0.25 


2 


09 


31 


55.2 


21 


28 03.3 


39.3 


5.4 


n A a 
U.44 


3 


09 


31 


57.7 


21 


30 16.5 


11.1 


2.5 


0.08 


4 


09 


31 


58.4 


21 


25 38.0 


8.8 


2.7 


0.10 


5 


09 


31 


59.0 


21 


24 25.3 


51.7 


6.2 


0.58 


6 


09 


31 


59.9 


21 


27 26.6 


51.8 


6.1 


0.58 


7 


09 


32 


00.6 


21 


23 38.5 


6.6 


2.5 


0.07 


8 


09 


32 


01.1 


21 


31 59.2 


8.7 


2.4 


0.06 


9 


09 


32 


01.1 


21 


32 33.7 


35.0 


5.0 


0.39 


10 


09 


32 


01.1 


21 


26 34.4 


12.4 


3.1 


0.09 


11 


09 


32 


01.2 


21 


30 46.5 


y.4 


O A 

2.4 


O.U 1 


12 


09 


32 


01.2 


21 


32 52.0 


189.1 


12.3 


1.40 


13 


09 


32 


01.9 


21 


31 11.1 


290.4 


15.7 


2.12 


14 


09 


32 


02.5 


21 


33 37.5 


9.4 


3.1 


0.07 


15 


09 


32 


03.0 


21 


26 55.1 


36.9 


5.3 


0.41 


16 


09 


32 


03.4 


21 


33 19.8 


18.5 


3.4 


0.14 


17 


09 


32 


03.7 


21 


29 32.9 


13.3 


3.2 


0.10 


18 


09 


32 


03.9 


21 


29 16.7 


21.1 


3.9 


0.15 


19 


09 


32 


04.6 


21 


29 59.9 


11.5 


3.0 


0.09 


20 


09 


32 


04.6 


21 


33 26.6 


318.0 


16.1 


2.35 


21 


09 


32 


05.0 


21 


26 41.4 


12.2 


2.5 


0.14 


22 


09 


32 


05.4 


21 


32 35.0 


1043.0 


29.5 


7.68 


23 


09 


32 


06.0 


21 


28 12.8 


157.1 


10.8 


1.13 


24 


09 


32 


06.1 


21 


28 18.3 


7.3 


2.5 


0.05 


25 


09 


32 


06.2 


21 


30 58.7 


2848.3 


48.4 


20.94 


26 


09 


32 


06.3 


21 


29 57.4 


44.5 


5.6 


0.33 


27 


09 


32 


06.6 


21 


30 04.8 


15.9 


3.5 


0.12 


28 


09 


32 


06.9 


21 


28 34.7 


9.9 


2.7 


0.07 


29 


09 


32 


06.9 


21 


28 55.0 


17.0 


2.9 


0.12 


30 


09 


32 


07.2 


21 


29 36.8 


47.9 


5.8 


0.34 


31 


09 


32 


07.4 


21 


30 55.1 


90.5 


8.5 


0.66 


32 


09 


32 


07.6 


21 


29 02.1 


40.7 


5.3 


0.29 


33 


09 


32 


07.7 


21 


29 56.2 


45.1 


5.6 


0.33 


34 


09 


32 


07.9 


21 


29 30.5 


30.6 


4.3 


0.22 


35 


09 


32 


07.9 


21 


31 44.9 


5.9 


2.6 


0.04 


36 


09 


32 


07.9 


21 


32 13.2 


31.2 


5.0 


0.23 


37 


09 


32 


09.1 


21 


29 09.0 


190.9 


11.9 


1.37 


38 


09 


32 


09.5 


21 


30 06.5 


Z4Z.U 


13.4 


1.79 


39 


09 


32 


09.6 


21 


28 49.1 


26.0 


A A 

4.4 


0.19 


40 


09 


32 


09.7 


21 


31 06.8 


1234.9 


32.0 


8.95 


41 


09 


32 


09.7 


21 


28 36.9 


28.4 


4.7 


0.20 


42 


09 


32 


09.7 


21 


24 24.1 


18.8 


3.6 


0.21 


43 


09 


32 


09.7 


21 


29 06.3 


152.1 


10.2 


1.09 


44 


09 


32 


09.8 


21 


30 05.8 


208.6 


10.6 


1.54 


45 


09 


32 


09.8 


21 


32 17.2 


47.0 


6.3 


0.35 


46 


09 


32 


09.9 


21 


34 32.5 


30.7 


A A 

4.4 


0.23 


47 


09 


32 


10.0 


21 


30 05.7 


140.3 


6.1 


1.04 


48 


09 


32 


10.0 


21 


29 52.0 


8.7 


2. 1 


0.06 


A n 
49 


09 


32 


10.1 


21 


30 01.2 


275.5 


13.1 


2.03 


50 


09 


32 


10.1 


21 


30 09.8 


103.0 


5.8 


1.15 


51 


09 


32 


10.1 


21 


30 08.5 


782.4 


22.4 


5.82 


52 


09 


32 


10.2 


21 


30 02.7 


93.7 


5.5 


0.69 


53 


09 


32 


10.4 


21 


29 58.6 


108.8 


8.5 


0.80 


54 


09 


32 


10.5 


21 


31 12.0 


9.6 


2.6 


0.07 


55 


09 


32 


10.8 


21 


31 26.4 


176.1 


12.0 


1.31 


56 


09 


32 


10.9 


21 


29 37.7 


12.1 


2.7 


0.09 


57 


09 


32 


11.5 


21 


27 41.7 


11.4 


3.1 


0.08 


58 


09 


32 


11.5 


21 


30 26.1 


137.9 


10.5 


1.02 


59 


09 


32 


11.8 


21 


28 34.9 


67.7 


/ .0 


0.49 


60 


09 


32 


11.8 


21 


24 57.3 


46.1 


5.6 


0.52 


61 


09 


32 


11.9 


21 


27 11.2 


29.5 


4.7 


0.21 


62 


09 


32 


12.0 


21 


30 33.9 


21.3 


4.3 


0.15 


63 


09 


32 


12.2 


21 


27 49.6 


40.0 


5.4 


0.29 


64 


09 


32 


12.3 


21 


29 22.9 


627.1 


21.6 


4.59 


65 


09 


32 


12.4 


21 


30 49.5 


11U.4 


9.6 


0.8i 


66 


09 


32 


12.5 


21 


30 25.9 


13.1 


3.1 


0.09 


67 


09 


32 


12.6 


21 


24 39.9 


346.9 


15.9 


3.91 


68 


09 


32 


12.6 


21 


32 08.1 


42.8 


6.1 


0.31 


69 


09 


32 


12.7 


21 


29 38.8 


8.8 


2. 1 


0.07 


70 


09 


32 


12.7 


21 


29 50.3 


34.7 


5.1 


0.26 


71 


09 


32 


13.1 


21 


30 55.8 


45.7 


6.1 


0.33 


72 


09 


32 


13.2 


21 


30 37.4 


31.1 


4.9 


0.23 


73 


09 


32 


13.6 


21 


35 33.0 


82.6 


7.2 


0.61 


74 


09 


32 


13.9 


21 


31 01.2 


17.1 


3.6 


0.13 


75 


09 


32 


14.0 


21 


33 05.2 


73.7 


7.6 


0.54 


76 


09 


32 


14.6 


21 


30 28.6 


40.3 


5.7 


0.30 


77 


09 


32 


14.6 


21 


31 01.2 


12.2 


3.3 


0.09 


78 


09 


32 


15.3 


21 


30 59.0 


5.3 


2.4 


0.04 


79 


09 


32 


15.4 


21 


29 24.9 


29.0 


4.5 


0.21 


80 


09 


32 


15.4 


21 


31 36.9 


15.2 


3.7 


0.11 
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TABLE 3 
Properties of X-Ray Point Sources 





R.A. 
(J2000) 


Decl. 
(J2000) 


Counts 


S/N 


L^/IO 38 
(erg s _1 ) 


81 


09 32 16.4 


21 31 11.7 


8.9 


2.4 


0.07 


82 


09 32 17.3 


21 29 46.5 


77.8 


7.9 


0.57 


83 


09 32 17.7 


21 30 20.4 


40.5 


5.6 


0.30 


84 


09 32 18.9 


21 26 36.1 


6.5 


2.4 


0.05 


85 


09 32 19.6 


21 31 05.3 


86.0 


8.3 


0.63 


86 


09 32 20.0 


21 27 48.9 


18.1 


3.8 


0.13 


87 


09 32 21.0 


21 30 02.2 


34.7 


5.5 


0.26 


88 


09 32 21.6 


21 28 20.9 


19.3 


4.0 


0.14 


89 


09 32 22.2 


21 32 41.8 


21.9 


3.7 


0.16 


90 


09 32 22.9 


21 29 24.7 


34.6 


5.2 


0.25 


91 


09 32 23.0 


21 33 02.5 


25.1 


4.1 


0.19 


92 


09 32 23.6 


21 29 13.7 


12.5 


2.9 


0.09 


a Intrinsic 
of index T 


luminosity (0.5-8.0 keV) 
= 1.95, n H =2.9x 10 20 cm" 


assuming a 

-2 


power law spectrum 



hot gas observed along lines of sight through the west 
zone is likely either (also) from the disk or from the wind 
located above the disk resulting in a relatively lower ab- 
sorption column density. This geometry naturally also 
explains the higher X-ray surface brightness to the west 
(Figure 2) but does not easily explain the observed hot 
gas morphology; the shape of the soft diffuse X-ray emis- 
sion does not approximate the bi-conical outflow seen in 
M82, for instance. 

In any case, the X-ray CCD data alone do not give ac- 
curate kinematic information of hot gas. The kinematics 
of an outflowing gas can be obtained from high- resolution 
spectroscopic measurements of warm ionized gas. In fact, 
Hagele et al. (2009) argue that H/3 line structures and 
properties measured from some nuclear star clusters are 
akin to those of Ha emission in NGC 1569 investigated 
by Westmoquette et al. (2007b, a) who concluded that 
these features are indicative of an outflow. Mapping the 
warm ionized gas in the central region of NGC 2903 will 
constrain whether or not the hot gas is escaping from the 
galaxy. 

5.3. Nuclear Activity 

There is no compact source detected coincident with 
the position of the mass center in NGC 2903. This 
is consistent with the lack of any AGN signature visi- 
ble at other wavelengths. We place an upper limit of 
Lx(o. 5-8.0 keV) = (1 — 2) x 10 38 erg s _1 to any undetected 
point-like nuclear source. 

We can infer the mass of the central object from the 
M B n - cr relation (Ferrarese & Merritt 2000). The ob- 
served stellar velocity dispersion in the nuclear region of 
NGC 2903 listed in the HYPERLEDA 3 database (Pa- 
turel et al. 2003) is 101.1±6.8 km s" 1 (Heraudeau & 
Simien 1998). This implies a black hole mass of 10 6 - 
10 7 M Q (Tremaine et al. 2002; Graham et al. 2011; Gra- 
ham 2012). Therefore, we take 10 7 M as an upper limit 
to the mass of the nucleus. If we assume that the am- 
bient hot gas is spherically accreting onto this central 
object, the volume-weighted Bondi accretion rate is es- 
timated as Meondi = 5xl0~ 8 M Q yr _1 using Equation 
(6) in Soria et al. (2006) and values of the gas tempera- 
ture and density as listed in Table 2. We emphasize that 

3 http://leda.univ-lyonl.fr 



this accretion rate only accounts for accretion of the hot 
phase of gas in the nuclear region. 

We define the dimensionless accretion parameter, m, as 
m = 0.1Mc 2 /iEdd = M/M^dd, where M is the accretion 
rate and Msad is the accretion rate that would produce 
the Eddington luminosity in a radiatively efficient case. 
For the NGC 2903 nucleus, we obtain rh = 2xl0~ 7 . In 
such low radiative efficiency regime, accretion is likely 
to be advection dominated. From the self-similar solu- 
tion of Narayan & Yi (1994), the bolometric luminosity 

is estimated as L DO i = f]Mc 2 = (10to)Mc 2 , where r] is 
the radiative efficiency. This leads to the corresponding 
Lbo\ ~10 34 erg s -1 . The predicted X-ray luminosity can 
be a factor of 10 lower than the bolometric luminosity for 
a supermassive black hole (Elvis et al. 1994), resulting in 
Lx ^10 33 erg s _1 . A similar Lx value is also obtained 
from the Lx-rh relation derived by Merloni et al. (2003) 
for advection models (taking into account that the al- 
ternative definition of m = Mc 2 / L^dd is used in their 
paper.) The hot gas component alone can only provide 
enough accretion to produce Lx ~10 33 erg s^ 1 . 

In general, nuclear black holes in nuclear starburst 
galaxies are much more active than in NGC 2903 be- 
cause cold gas is also accreted onto the black hole. There 
is a tight correlation between nuclear SFR and black 
hole accretion rate (BHAR) found in Seyfert galaxies 
(Diamond-Stanic & Rieke 2012). The SFR-BHAR sug- 
gests that external cold gas triggers star formation and 
AGNs or nuclear star formation fuels subsequent AGN 
activity. The SFR measured from hydrogen recombina- 
tion lines in the central 15" (650 pc) region in NGC 2903 
is 0.7 M & yr _1 (Alonso-Herrero et al. 2001). For such 
an SFR, we would expect M ~ 0.01-0.1 M Q yr" 1 for 
the Seyfert galaxies studied by (Diamond-Stanic & Rieke 
2012) because the SFR corresponds to a cold gas inflow 
M. However, given the upper limit to the X-ray lu- 
minosity, the cold gas M in NGC 2903 must be much 
lower, < 3xl0 -5 M Q yr _1 , and is inconsistent with the 
SFR-BHAR correlation, even considering systematic un- 
certainties among different methods to derive SFRs and 
accretion rates. 

Planesas et al. (1997) estimated an H2 mass of 
1.8x10 s Mq, corresponding to a surface density of 
160 M Q pc" 2 in the nuclear region of NGC 2903. More- 
over, the recent high resolution CO image shows a gas 
structure swirling into the center (P.-Y. Hsieh, private 
communication). Therefore, there is enough material 
available in the central region to accrete on to a central 
massive object in NGC 2903. 

We point out that only a small fraction (~ 10%) of 
gas fueling nuclear SFR is needed in order to trigger a 
nuclear activity as luminous as in Seyfert galaxies. In 
NGC 2903, the 10% of cold gas inflow may be also con- 
sumed in the nuclear starburst, heated by massive stars 
and supernovae, and/or ejected in a galactic outflow. 



Support for this work was provided in part by the 
National Aeronautics and Space Administration through 
Chandra Award GO0-11099A issued by the Chandra X- 
ray Observatory Center, which is operated by the Smith- 
sonian Astrophysical Observatory for and on behalf of 
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TABLE 4 

Fit Parameters of Bright X-Ray Point Sources 



Parameter 


J09325.4+213235 


J093206. 2+213058 


J093209. 7+213106 


J093210.1+213008 


J093212.3+212922 




McKaL 


n H /10 22 (cm- 2 ) 
kT (keV) 

z/z e 

K/10-* 
x 2 /dof 

/xo.5-8.0kev/10 14 (erg s 1 cm 2 ) 


15+ - 02 
u - lo -0.02 

o.uo_ 2g 

0.82lHr 
88.9/77 
6.67 


14+ 001 

3 n 4 +0.24 
J - u -0.24 
05+ 011 
o i 7+0.14 

148/156 
18.2 


n 90+0.03 

u - zo -0.03 

1 94+0.19 

^ -0.17 
03+ ' 05 

, ocr+0.14 
L - io -0.14 

64.1/90 
7.04 


2 00+ ' 24 
z.uu_ 25 

9 O1+0.35 
z -°-0.34 

14+ 025 
u -l^-0.14 

9 9Q+O.O7 

68.6/67 
7.44 


12+ - 04 
u - ±z -0.04 
7 01 +3- 15 
' •° 1 -1.87 
00+ ' 57 

0.38t°'f 
81.5/50 
4.74 


Power law 


n H /10 22 (cm" 2 ) 

r 

K/10- 4 
X 2 /dof 

/xo.5-8.0kev/10 14 (erg s 1 cm 2 ) 


25+ - 02 

2 26+ 007 

27+ 03 
u - z ' -0.03 

101/78 
6.82 


25+ 02 
u - zo -0.02 

2 29+ - 07 
z - zy -0.06 
7C :+0.05 
u - '°-0.05 

151/157 
18.8 


40+ ' 03 
u - 4u -0.02 
9 09+0. 09 
z -° -0.09 

73.3/91 
7.17 


4+0.29 

1 RR+ ' 87 
1-O»_0 .89 

70.9/68 
7.59 


n , „+0.03 

0.16_ 004 

1-6418:81 
o.o9i8:8i 

85.6/51 
5.00 


DiskBB 


n H /10 22 (cm- 2 ) 
kT in (keV) 
AT/10- 4 
X 2 /dof 

/xo.5-8.0kev/10 14 (erg s 1 cm 2 ) 


ri5+ - 02 

U - UO -0.02 

i-07tg:°0B 7 
28-i±S:8 

84.5/78 
6.59 


U - UO -0.01 

1 oi+ - 04 

91 4+ 25 ' 1 
<"■■*— 12.9 

171/157 
17.3 


0.12t° ;°| 

04 0+24.4 

68.9/91 
6.89 


1 61+ ' 38 

!- D1 _0.32 

1 ns+ - 17 

i - Ucl -0.14 

67 5+ 37 ' 8 
D ' -°-23.0 

68.2/68 
7.00 


04+ ' 03 

1 41+0-17 

e.osl 2 -^ 

73.5/51 
4.42 


X-ray luminosity 


^xrj.s-s.okev/ 1038 ( cr § s_1 ) 


6.64 


20.8 


9.50 


13.8 


4.33 



Note. — Errors are lcr (68% confidence interval) uncertainty. Intrinsic X-ray luminosity is estimated using the best-fit model. 
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APPENDIX 
X-RAY POINT-SOURCE ANALYSIS 

Source positions, number of source counts (corrected for the finite-aperture model PSF), S/N, and estimated source 
luminosities in the 0.5—8.0 keV range are tabulated in Table 3. The listed source positions are refinements from 
the initial source detection estimates; refined positions were made by fitting circular Gaussian models to the spatial 
distribution of X-ray events in the vicinity of each source. The source luminosities were estimated from the average 
count rate during the observation using the Portable Interactive Multi-Mission Simulator (PIMMS) 4 assuming an 
absorbed power-law spectral shape with spectral index T = 1.95 and with a hydrogen column density equal to the 
Galactic column density along the line of sight of NGC 2903, nn = 2.9 x 10 20 cm -2 , and solar abundances. We 
note that estimating flux using individual response and effective area files for each source gives more accurate value. 
However, PIMMS estimation is reasonable here due to low S/N for the most of sources; therefore, the statistical 
uncertainties are larger than the systematic uncertainties due to differences in the response matrices implemented in 
PIMMS than those of the actual observation. 

Here, the X-ray spectra of the five highest-count sources are analyzed. This includes one source in the nuclear 
region, CXOUJ093210. 1+213008, the source mentioned in § 2.1, and three others. All five are coincident with sources 
detected with XMM-Newton in a 21 April 2009 observation (Perez-Ramirez et al. 2010) although sources in the 
nuclear region are highly confused at the XMM-Newton resolution (XMM-NGC2903 X-l is 3.9" distant from CX- 
OUJ093210.1+213008). 

Spectral analysis was performed using the XSPEC v. 12. 6 spectral-fitting package using redistribution matrices and 
ancillary response files appropriate to the source location on the detector. These files were generated using the CIAO 
utilities mkacisrmf, mkarf, and associated programs. Absorbed (xSPEC's phabs) power law (powerlaw), blackbody 
accretion disk (diskbb), and optically-thin thermal plasma (mekal) models were fitted to the observed spectra in the 
0.3—10.0 keV range using the x 2 fit statistic (events were grouped to ensure a minimum of 10 counts per spectral 
energy bin) . Results of the fitting are given in Table 4 and the observed and model spectra and delta \ 2 are displayed 
in Figure 9. The intrinsic (absorption-corrected) 0.5—8.0 keV luminosities of these sources, as estimated using the 
best-fitting model, are listed in the bottom row of Table 4. 

4 http://heasarc.nasa.gov/docs/softwarc/tools/pimms.html 
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Three sources are best fit using the blackbody-disk model. This model is typically applied to accreting X-ray binary 
systems. The diskbb model results in inner disk temperatures of T; n ~ 1 keV as expected for stellar-mass black holes 
(e.g., Remillard & McClintock 2006). The most luminous source, CXOUJ093206. 2+213058, and the source located in 
the bright star-forming region, CXOUJ093209. 7+213106, are best fit using a thermal emission model indicative of hot 
gas. Note, however, that the best-fitting abundances for these (and for the other three) sources are much less than 
the solar value due to the lack of observed line features. Furthermore, the resulting plasma temperatures are much 
higher than is typical of diffuse hot gas in nearby spiral galaxies such as NGC 2903. These results suggest that these 
two sources may also be accreting X-ray binaries. The source in the nuclear region has by far the highest intervening 
column density. This absorbing material is likely local to the nuclear region. 
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Fig. 9. — Spectra, best-fitting models (solid histogram), and fit residuals (lower panel) of the brightest sources in NGC 2903 detected in 
the Chandra observation taken on March 7, 20f 0. Details of all model fits including fit parameters and X-ray fluxes estimated from model 
fits arc listed in Table 4. 



